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The  transformations  of  Gaussian  radiation  beams  caused  by  reflection 
off  mirrors  is  an  important  issue  for  FELs  operating  as  oscillators.  The 
reflected  radiation  from  a  single  incident  Gaussian  mode  will  contain  other 
modes  due  to  the  finite  mirror  size,  the  deflection  of  the  beam  and 
mismatches  in  the  curvature.  A  method  for  analytic  computation  of  the 
reflection  matrix  is  developed  by  taking  the  convolution  of  the  source 
function  at  the  surface  of  the  mirror  with  the  paraxial  propagator.  The 
mirror  surface  that  reflects  spherical  incoming  wavefronts  into  spherical 
outgoing  is  found  to  be  a  paraboloid.  Integral^  expressions  for  the 
V  reflection  coef  f icients  (r^J  \  for  any  incoming  mode  u^nX  iid:o  the  outgoing 

are  obtained^as  functions  '"of  the  deflection  angle  the  reflected  beam 
spot  size  V  and  the  mirror  size.  The  coefficient  R^~ 


igle  (tf  th 


for  the  lovest-to- 


lowest  mode  reflection  is  determined  analytically. 


then  be  selected,  depending  on  the  mirror  size,  to  maximizeVRr 


00,  V 

The  spot  size  W-  can 
rva  '°7 


The  ratio 


of  the  mirror  size  to  the  spot  size  is  the  dominant  factor  determining  the 
reflection  coefficient.  The  effects  of  deflecting  the  light  beam  enter  as 
small  corrections,  of  first  order  in  the  diffraction  angle,  9^  «  1.  _ 
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I.  INTRODUCTION 


Free  Electron  Lasers  (PEL)  operating  as  oscillators*-^  require  the 

8-10 

trapping  of  light  pulses  between  systems  of  mirrors  (resonators). 

These  pulses  are  repeatedly  amplified  via  synchronous  interaction  with 

electron  pulses  passing  through  the  wiggler.  The  simplest  oscillator 

configuration  is  that  of  an  open  resonator  with  two  opposed  identical 

mirrors.  The  radiation  vector  potential  for  this  arrangement  is  expressed 

in  terms  of  the  free  space  eigenmodes  A  (r)  =  A  (r)e  of  the  paraxial 

an  mn  an 

wave  equation**,  where  is  the  polarization  vector  and 


A(D  = 

mn 


Umn(r;W)  i|k2+ 

2. 1/2  e  L 


2R(z) 


1Smn<2> 


**h) 


The  exponent  *(r)  *  k[z  +  (x  +y^)/2R(z) ]  in  (1)  contains  the  rapidly 


varying  phase  on  the  wavelength  scale  X  =  2n/k. 


The  constant  phase 


2  2 

wavefronts  are  spherical  of  curvature  1/R(z)  =  z/(z  +  b  ).  The  spot  size 

?  7  1/7  1/7 

of  the  radiation  envelope  is  W(z)  =  w  (1  +  z  /b  )  ,  where  w  =  (2b/k) 

is  the  waist,  and  the  distance  z  is  measured  from  the  position  of  the 
waist.  The  amplitude  squared  of  the  mode  drops  by  1/2  over  a  distance 
equal  to  the  Rayleigh  length  b  (also  known  as  confocal  parameter).  Most  of 
the  radiation  is  confined  within  a  cone  parametrized  by  the  diffraction 


angle  9^  =  W/z  =  (X/bn)*"‘‘.  The  amplitude  profile  umn(r;V)  contains  the 
transverse  spatial  variation,  equivalent  to  a  small  k^,  perpendicular  to 
the  z-direction.  Higher  modes  correspond  to  a  larger  effective  k^, 
responsible  for  the  phase  shift  exp  (i8mn(z)].  For  a  given  k  the  mode  is 
completely  defined  by  the  two  independent  parameters  R  and  w  (or  any 


combination  of  two  out  of  the  four  quantities  R,  w,  z  and  b). 


The  functions  u|nn(r;W)  depend  on  the  elected  coordinate  system.  In 
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rectangular  coordinates  (x,y,z)  they  are  given  by 


(t)  "n^)  *  '  ^  ■  <2> 

where  H  ,  H  are  the  Hermite  polynomials,  a  =  (2/V^)^/'^(n  2m+n  m!  n!)_1/^ 
m  n  mn 

is  the  normalization  factor,  setting  the  total  electromagnetic  flux  carried 
by  the  mode  equal  to  unity,  and  Smn(z)  =  (m  +  n  +  1)  tan~*(z/b)  is  the  slow 
phase.  In  cylindrical  coordinates  (r,0,z),  where  tan©  =  x/y,  r  = 
(x^+y^)*^,  u^(r,0;V)  take  the  form 


(3) 


where  +p  (-p)  signifies  cosine  (sine)  poloidal  dependence,  are  the 
Laguerre  polynomials,  a^  =  (l/2ltW^)^^  [m!/(m  +  p)!] and  6^(z)  =  (2m  + 

p  +  1)  tan~*(z/b). 

The  electron  beam  is  an  optically  active  medium  that  alters  the 

characteristic  parameters  of  the  radiation  after  each  passage.  During  the 

build-up  period  the  modal  content  and  the  structure  of  the  light  pulses 

inside  the  oscillator  will  change.  A  numerical  method  has  been  developed 

recently  optimizing  the  representation  for  the  amplified  radiation.  In  the 
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source  dependent  expansion  ’  the  waist  size  and  the  curvature  of  the 
elected  modal  basis  is  tailored  according  to  the  driving  source  term.  That 
minimizes  the  number  of  modes  required  to  describe  the  light  beam.  In 
general,  the  curvature  and  waist  size  of  these  modes  does  not  match  the 
curvature  and  waist  of  the  vacuum  eigenmodes  for  the  resonator.  Therefore 
the  transfer  matrix  for  a  given  mirror  must  be  known  for  arbitrary  incoming 
modes.  This  need  stems  from  computational  as  well  as  physical  reasons. 
The  knowledge  of  the  cavity  reflection  matrix  R,  together  with  the  gain 


2 


matrix  G  through  the  viggler,  is  necessary  in  determining  the  potential  for 
steady  state  operation. 

During  high  power  operation,  grazing  mirror  incidence  may  be  necessary 
to  avoid  exceeding  the  dielectric  breakdown  limit  for  the  reflection 
surface.  Also,  in  case  of  a  high  per-pass  gain  with  optical  guiding,  the 
waist  for  the  reflected  radiation  would  be  much  larger  than  the  incoming. 
In  two-mirror  resonators  the  reflected  radiation  could  hit  the  wiggler. 
Therefore,  ring  resonators,  including  three  or  more  mirrors,  must  be 
employed  for  the  deflection  and  recirculation  of  the  light  pulses.  The 
study  of  the  reflection  matrix  must  be  extended  to  include  the  effects  of 
deflecting  the  light  beam,  in  addition  to  finite  mirror  size  and  curvature 
mismatches. 

The  amplitude  profile  of  the  incoming  radiation  will  be  modified  by 
reflection.  A  single  incident  mode  Amn(r^)  will,  in  general,  be  partially 
reflected  into  different  modes  Apq(rQ)  where  (m,n)  *  (p,q).  This  is  caused 
by  the  deflection  of  the  light  beam,  the  finite  size  of  the  mirror  and  the 
curvature  mismatches.  Reflection  into  other  modes  will  affect  the 
interaction  between  the  electron  beam  and  the  radiation  in  a  number  of 
ways.  First  the  rms  radius  of  the  light  beam  will  change,  affecting  the 
matching  beam  condition.  Second,  the  light  pulse  will  spread  axially 
because  of  dispersion  among  different  modes,  since  the  phase  velocity 
depends  on  the  modal  number  (m,n).  Third,  different  phase  advances  during 
reflection  among  the  various  modes  may  render  these  modes  out  of  phase 
after  a  number  of  bouncings  off  the  resonator.  7or  the  above  reasons  the 
fraction  of  radiation  scattered  into  other  modes  will  contribute  to  the 
losses  in  FEL  oscillators. 
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II.  OUTLINE  OF  THE  METHOD 


The  geometry  of  the  reflection  is  illustrated  in  Fig.  1.  The 

subscripts  i  and  o  denote  the  coordinate  system  used  for  expressing 

incoming  and  outgoing  modes.  r^  is  defined  with  the  z^  axis  along  the 

direction  of  incidence  and  r  has  the  z  axis  along  the  direction  of 

o  o 

reflection.  The  origins  are  displaced  from  the  mirror  center  by  lj  and  1 

respectively,  where  1^  is  the  distance  of  the  minimum  waist  w^  for  the 

incoming  radiation  and  1  is  the  distance  of  the  minimum  waist  w  for  the 

o  o 

reflected  mode.  A  third  coordinate  system  rg  with  the  origin  at  the  mirror 
center  and  zg  aligned  with  the  mirror  axis  will  be  useful  in  the 
computations.  Underlined  quantities  r^,  rQ  and  £s  stand  for  the  mirror 
surface  coordinates  in  each  reference  frame.  The  transformations  among  the 
various  frames  are  defined  by 


X.  =  x  cos  £  - 
is  l 


Z  Sln  O 
s  ^ 


yi  =  ys’ 


<4a>  y0  =  ys’ 


(4b) 


z . 
1 


z  cos 
s 


2  +  xs 


sin  1+1. 


z 

o 


z  cos  ^  -  x  sin  ^  +  1  . 
s  z  s  z  o 


We  are  interested  in  cases  when  the  reflected  radiation  remains 

focused  along  some  direction  zq  making  an  angle  $  with  the  incoming  z^. 

Then  the  reflected  vector  potential  will  also  be  expandable  in  free  space 

eigenmodes  Ap^(ro>  of  the  paraxial  equation  in  the  new  direction.  The 

mirror  surface  generating  focused  reflection  in  the  desired  direction  can 

not  be  arbitrary  but  must  be  appropriately  defined.  The  angle  of 

deflection  $  will  enter  the  equation  defining  the  mirror  surface.  The  other 

surface  parameter,  namely  the  curvature  1/R  ,  is  a  free  parameter.  It 

m 

determines  the  curvature  1/Rq  for  the  outgoing  modes  given  the  curvature 
1/Rj  of  the  incoming  modes.  In  case  of  reflection  by  an  arbitrary  mirror 
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surface,  the  scattered  radiation  cannot,  in  general,  be  covered  by  the 
paraxial  modes  that  do  not  form  a  complete  set  in  three  dimensions. 

We  consider  incoming  radiation  of  given  curvature  and  of  arbitrary 
amplitude  profile  A1(rj),  consisting  of  various  modes  (m,n)  vith  the  same 
R^(z).  If  both  incident  and  reflected  radiation  are  expanded  into 
eigenmodes, 


A1(r.)  =  V  c1  A  (r.)  , 

'  i '  u  mn  mn'  r 


m,n 


-  If  -«  W  • 


(5) 


the  relation  among  the  incident  and  reflected  expansion  coefficients  c1^, 

c°  is  formulated  by 

pq  f 

c°=  R  c*  ,  (6a) 


or 


pq 


=  E  R1 


m,n 


mn 

pq 


mn 


(6b) 


where  Rmn  are  the  elements  of  the  reflection  matrix  R. 

pq 

We  examine  the  case  when  the  mirror  dimensions  p  are  much  larger  than 

the  wavelength  A,  A  «  p  (otherwise  diffraction  rather  than  reflection 

would  prevail).  We  also  assume  that  the  angle  £  subtended  by  the  mirror 

C  =  p/R  ,  where  R  parametrizes  the  radius  of  curvature,  is  small,  of  the 
m  m 

order  of  the  diffraction  angle  0^,  C  ~  0^  ~  e.  The  j-th  component  of  the 
reflected  vector  potential  at  distance  |rQ  -  rQ|  »  A  from  the  mirror 
surface  S  is  then  given  by 


A°  (r  )  _  _  — 

A  (v)vV  “  2 n 


ik  r  -  r 
e  1  o  -o' 


ds 


'V  Eol 


(v)'-o 


(r  )  (n  •  Ar) .  (7) 
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In  Eq.  (7)  n*Ar  is  the  obliqueness  factor  where  Ar  =  (r  -r  )/|r  -r  I  and  n 

o  — o  1  o  *-0 1 

is  the  normal  unit  vector  to  the  reflecting  surface.  The  surface  element 
ds  is  given  by  ds  =  Slzo-f(xo,yo)]dxQdyodzo  where  zq  =  f(xQ,yo)  is  the 
surface  equation.  Equation  (7)  is  the  convolution  of  a  source  term  As(rQ) 
at  the  mirror  surface  with  the  propagator  exp(ik lr0-rQ | )/ |rQ-ro | ,  i.e.,  a 
superposition  of  spherical  waves  originating  at  S.  The  source  term  As(rQ) 
is  specified  from  the  incoming  vector  potential  A1(rj)  through  the  boundary 
conditions  and  the  coordinate  transformations  (4).  We  will  assume  a 
perfectly  conducting  surface,  where  the  incident  and  reflected  fields  are 
related  by 

As  =  -A1  +  2  (n  •  A1)  n  ,  (8a) 


and  n  is  the  normal  unit  vector  to  the  reflecting  surface.  When  the  solid 
angle  subtended  by  the  mirror  is  small,  the  last  term  in  (8a)  is  very  small 
and  the  boundary  condition  becomes, 


A 


s 


(V) 


(8b) 


Relation  (8b)  corresponds  to  a  phase  shift  by  n  during  reflection.  It  is 

independent  of  the  wave  polarization,  thus  the  subscript  (\>)  is  dropped. 

Cross  polarization  effects,  due  to  the  last  term  in  (8a)  are  discussed  in 

Ref.  10.  Most  of  the  computations  will  be  performed  on  the  mirror  surface. 

To  simplify  the  notation  from  now  on,  we  drop  the  bar  (_)  under  the  mirror 

coordinates  r.  Subscripted  quantities  such  as  r^,  rQ,  rg  will  signify  the 

mirror  surface  in  each  reference  frame.  Unsubscripted  coordinates  will 

denote  the  observation  point  in  the  reflected  radiation  frame  of  reference. 

We  seek  cases  when  the  reflected  radiation  propagates  focused  along 

z  .  contained  within  a  crdss  section  of  dimensions  x,y  «  z-z  .  The 
o  o 
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expansion  |r-rQ|  -  (z-zq)  {1  +  [(x-xq)  +  (y-yQ)  ]/2(z-zq)  }  replaces  the 
full  propagator  inside  (7)  with  the  paraxial  propagator  in  that 
direction, 


-ik(z-zQ) 


(r,  rQ)  = 


<x-xq)2+  (y-yQ)2 

2<z"zo) 


It  is  known  that  the  profile  of  a  given  eigenmode  Amn(xo,yQ,zo)  at  zq 

is  generated  by  the  propagator  Uj^^r  )  acting  on  the  mode  Amn(x,y,0)  at 

z=0.  The  inverse  propagator  U^(r,ro)  therefore  reproduces  Amn(x,y,0)  from 

Amn^xo,yo,zo^ '  Th*s  suggests  multiplying  and  dividing  the  integrand  inside 

j  9  1/2 

(7)  by  exp  [i*(rQ)]  /  [1  +  z^/b^  ) 4 ,  recasting  (7)  in  the  form, 


iA(r  ) 


rr  v 

A°(r)  =  Jjds  e  °  S(rQ)  e 


i#  (r  ) 

'  r*  ' 


U-k(r’ro>  - 


where  the  source  term  S(rQ)  is, 


z  2(r  )11/2 

S(rQ)  =  Ai[ri(ro)]  (n  •  Ar)  1  +  0  2g 

b 

o 


and  the  phase  A(rQ)  is  given,  in  outgoing  coordinates  rQ,  by 

2,  N  2  ,.22, 
xi  (ro)  +  yi  (ro)  xo  +  yo 

A(r  )  =  k  z.(r  )  +  z  +  - 2 - 1 - +  _2 - #  (12) 

L  1  °  2Ri(ro)  2Ro(ro) 

The  expression  for  A(rQ)  depends  on  the  angle  ♦  through  the  transformations 
between  the  incoming  and  the  outgoing  coordinates,  Eqs.  (4).  Expression 


(10)  is  the  approximation  of  the  exact  solution  (7)  to  order 
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I(x-xq)^  +  (y-yo)^]/2(z-zo)2  -  e^.  it  is  valid  provided  that  the  surface  S 
produces  focused  reflection  along  the  desired  direction.  Otherwise,  the 
paraxial  limit  will  fail  to  encompass  all  the  radiation  contained  in  the 
original  expression  (7). 

The  term  exp[id(rQ)J  is  varying  rapidly,  on  the  scale  of  the 
wavelength  X.  Therefore,  its  convolution  with  the  slowly  varying  source 
term  over  an  arbitrary  surface  will  be  vanishingly  small.  In  general,  this 
corresponds  to  radiation  scattering  where  only  a  small  fraction  of  the 
incoming  radiation  is  reflected  along  the  considered  direction  $.  The 
integral  (10)  will  be  finite  only  when  it  is  possible  to  satisfy  the 
condition  A(rQ)  =  constant  over  some  surface  S.  If,  in  addition,  S  is  much 
larger  than  X,  expression  (10)  will  be  finite  only  within  a  narrow  angle  6$ 
around  This  guarantees  that  the  reflected  radiation  remains  focused 
along  that  direction.  Therefore,  a  condition  that  the  exact  reflected 
radiation  (7)  be  fully  covered  by  the  paraxial  limit  (10)  is  that 

A(rQ)  =  const.,  (13) 


along  the  surface  S.  Accordingly,  the  optical  path  is  the  same  along  the 
various  rays  connecting  an  incoming  wavefront  with  its  mirror  image 
(reflected)  wavefront. 

Requirement  (13)  defines  the  appropriate  mirror  surface  zo=fQ(xo,yo; $) 

for  reflection  in  the  elected  direction.  Expressing  all  quantities  inside 

(13)  in  the  mirror  coordinate  frame,  applying  the  transformations  (4)  and 

2 

using  the  scaling  x  /R  -  y  /R  -  e  <<  1,  z  /R  -  e  we  obtain  from  (13) 
sm  sm  sm 


z 


s 


1 

2Rcos 

m 


1 

2 


cos 


2  * 


(14a) 


where 
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1 

R 


(14b) 


Equation  (14a)  is  the  analytic  expression  for  a  paraboloid  surface. 

Rm  parametrizes  the  mirror  curvature,  being  positive  or  negative  for  a 

convex  or  concave  mirror  respectively.  The  surface  is  reflection  symmetric 

with  (zx)  and  (zy)  as  the  symmetry  planes;  there  is  no  rotational 
s  s 

symmetry  around  zg.  Surface  (14a)  can  also  be  approximated,  to  second 
2  2 

order  in  (x  /R  )  ,  (y  /R  )  by  hyperboloids  or  ellipsoids  defined  by 
s  m  s  m 

(zs  -  Vos  !)  -  xs2cos2  2  -  ys2  *  Rm2cos2  2  ’  (15a) 

fz  +  R  cos  £)  +■  x  ^cos^  |  +  y  ^  =  R  ^cos^  f  .  (15b) 

v,  s  m  4/  s  z  s  m  z 


All  the  surfaces  become  spherical  in  the  limit  of  perpendicular  incidence  4> 
=  0,  and  plane  mirrors  when  Rm  -»  ®.  Relation  (14b)  defines  the  curvature 
of  the  reflected  modes  from  the  incoming  mode  curvature  and  the  curvature 
of  the  mirror. 

Switching  Eq.  (12)  into  the  mirror-aligned  coordinates  r  through  Eqs. 

s 

(4),  and  using  the  surface  constraints  (14),  it  follows  that 


P  >2. 


A(r  )  =  const, 
s 


°Mr)] 


where  p  parametrizes  the  mirror  size.  A  more  complicated  surface  equation 
(higher  than  quadratic  in  x,  y,  z)  is  required  to  improve  the  constancy  to 
a  higher  order.  Since  kp  »  1,  the  approximation  A(rg)  =  const,  is 
satisfactory  for  a  first  order  expansion  of  the  reflection  matrix  in  powers 
of  p/R  ,  as  long  as  p/R  <  (kp)-*.  In  case  p/R_  >  (kp)-*,  the  slow 
variation  of  A(xg,yg)  over  S  must  be  included.  That  introduces  an 
additional  contribution  in  the  reflection  matrix,  known  as  spherical 
aberation. 
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III.  COMPUTATION  OF  THE  REFLECTION  MATRIX 


The  reflected  radiation  is  expressed  by 


A°(r)  -  J  Jdx0dyo  ff(xo’yo)  e  000  U-k(r’ro) 

S 

where  a(xQ,yo)  *  S[xQ,yo,zo(xo,yo)] .  Expanding  the 
terms  of  umn(xo,yo), 


(16) 

source  a(x  ,y  )  in 


*(* 0»y0)  =  E  Rmn «. 

m ,  n 


(x  ,y  ;W  ). 
mn  o  o  o 


(17) 


and  exploiting  the  property  of  the  inverse  propagator  U  ^  the  reflected 
vector  potential  A°(r)  at  z  =  0  becomes 


A°(x,y,0)  =  £  Rmn  u  (x.yjW  ),  (18) 

m,n 

where  VQ(z)  =  w  (1  +  z^/b  ^)*^,  wq  =  (2b  /k)*^.  Expression  (18)  is  a 

complete  decomposition  of  the  reflected  radiation  into  paraxial  eigenmodes 

for  incident  radiation  of  arbitrary  profile. 

According  to  the  definition  (6b),  the  R™^  element  of  the  reflection 

matrix  R  is  obtained  from  the  source  term  ffpq(x0*yo)  inside  (16)  generated 

by  a  single  incident  eigenmode  Apq I r £ ( r0) I •  The  integration  is  performed 

in  the  mirror-aligned  coordinates,  taking  advantage  of  the  existing 

symmetries.  The  coordinates  and  r  ,  defining  the  incoming  and  outgoing 

wave  functions,  become  explicit  functions  of  xg,  yg  through  the 

transformations  (A).  The  surface  equation  (14a)  is  used  to  express  zg  in 

terms  of  (x  ,y  ).  The  mirror  boundary 
s  s 
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(19) 


2  2  * 
x„  cos  f  + 
s  l 


2 

P 


is  defined  by  the  intersection  of  the  infinite  surface  (14a)  with  the  plane 

2  2 

zg  =  const  =  2  p  cos  ($/2)/  R^.  After  the  above  manipulations  the 

reflection  matrix  elements  take  the  form 


^mn 

pq 


u  (x  ,y  )  u  (x  ,y  ) 
dx  dy  mn  5  5  pq  s--  s- 


s  '  s 


2  , 


1 


1/2 


-.1/2 


zo  (Vys> 


z.  (x  ,y  ) 
1  s  Js 


b.‘ 

1 


lA(xs,ys) 


(20) 


x  e 


iS1  (x  ,y  )  -  i6°  (x  ,y  ) 
pq  s  s  mn  s  Js 


2  .  2  *  , 

*  (.  xs  .  +  xs  sin  2  1 

:os  f  (1  -  r  sm  f  -  —2 - 


where 


u  (x  ,y  )  a  u  x  (x  ,y  ) ,  y  , 
mnv  s  s  mnL  o'  s,Js/’  •'sj  ’ 


u  (x  ,y  ) 
pq  s’-'s' 


pq 


2  2 

Expression  (20)  is  correct  to  order  p  /Rm  . 

Each  representation  of  R  is  tied  to  the  choice  of  the  basis  functions 
umn(r).  In  any  case  R,  as  given  by  (20),  depends  on  four  parameters^ 

R  *  R(*,  a,  u;  l).  (21) 


♦  is  the  reflection  angle  shown  in  Fig.  1.  a  is  the  ratio  of  the  incoming 

to  the  outgoing  spot  size  at  the  mirror,  a  =  W^(l.. )/¥q(1o) .  y  =  p/WQ 

parametrizes  the  mirror  size  compared  to  the  radiation  spot  size.  (,  = 

V  /R  scales  as  the  diffraction  angle  0,  =  W  /I  multiplied  by  the 

om  6  d  o  o  r  ' 

curvature  mismatch  R  /R  between  the  mirror  and  the  radiation  wavefronts. 

o  m 

The  spot  size  Wq  enters  as  a  free  parameter  because  only  the  curvature  1/Rq 
for  the  reflected  modes  is  specified  by  the  mirror  geometry.  Since  many 
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combinations  of  Wq  and  1  apply  to  a  given  curvature  according  to  paragraph 

Eq.  (1),  an  additional  selection  rule  for  is  needed.  Note  that  WQ  does 

not  have  to  match  W^.  This  is  obvious  in  cases  when  the  mirror  size  p  is 

smaller  that  V..  Each  value  of  V  defines  a  complete  set  of  modes  for  the 
1  o 

reflected  radiation  and  an  equivalent  representation  for  R. 

Parameters  ♦,  a,  and  p  can  be  arbitrary.  In  most  cases  of  interest, 
however,  £  is  small,  £  «  1,  of  the  same  order  as  the  diffraction  angle  0^. 
The  analytic  computation  of  the  matrix  elements  is  carried  out  by  expanding 
the  integral  (20)  in  powers  of  £, 

R  =  R(0)  +  l  R(l)  +  l2  R(2) .  (22) 

The  first  order  expansion  is  performed  in  Ref.  10.  In  this  paper  we 
review  some  of  the  general  properties  of  R  and  focus  on  the  reflection 
of  the  lowest  mode  Uqq. 


IV.  LIMITING  CASES 

When  the  mirror  radius  tends  to  infinity  (1/R^  *  0),  or  in  cases 
of  vertical  incidence  on  the  mirror  ($  =  0),  the  higher  order 
corrections  in  the  reflection  matrix  R  disappear, 

R  =  R(0)  (23) 

in  both  representations.  The  nondiagonal  elements  in  R  stem  from  the 
finite  mirror  size  only.  If,  in  addition,  the  mirror  size  is  very 
large,  p  »  1,  it  is  appropriate  to  take  VQ  =  as  best  representation 
for  the  reflected  radiation.  The  a  =  1  limit  yields 


(24) 


Thus,  in  case  of  large  curved  mirror  and  vertical  incidence,  or  large 
plane  mirror  and  arbitrary  incidence,  the  reflection  matrix  is  the 
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identity  matrix. 

The  case  a  =  1  is  of  special  interest  for  arbitrary  angle  of 

deflection  $  and  mirror  curvature  1/R,  as  it  will  be  explained  in  the 
next  section.  For  finite  mirror  size  p  >  VQ,  (  u  >1),  there  exists 
zeroth  order  nondiagonal  terms  inside  R(0).  Since  R(0)  is  independent 
of  the  angle  of  deflection  the  mirror  size  yields  the  dominant 

contribution  to  the  reflection  into  modes  different  than  the  incoming. 
The  effects  of  the  deflection  of  the  light  beam  enter  to  first  order  in 
l,  R(l) ,  or  higher. 

As  the  mirror  size  becomes  very  large  and  the  limits  of  integration 
are  extended  to  infinity  the  orthogonality  among  the  various  modes  u 

mn 

becomes  effective.  The  off-diagonal  terms  in  R(0)  become  comparable  to 

2 

the  first  order  terms  R(l)  roughly  when  1/p  ~  Z  ~  8^.  At  the  limit  p  -> 

”  all  the  nondiagonal  elements  of  R  are  reduced  to  order  £  or  higher, 

Rpq  =  *  Rpq<D  +  n  *P,  n  *  q,  (25a) 

and  the  only  matrix  elements  of  zeroth  order  in  ?,  are  the  diagonal 


Rmn 

mn 


Rmn 

mn 


(0)  +  oa2). 


(25b) 


in  both  Hermite  and  Laguerre  represenations.  The  lowest  correction  in 

2 

the  diagonal  elements  is  of  second  order  while  the  first  order 
contribution  disappears.  This  is  consistent  with  flux  conservation 
during  reflection  in  case  of  a  large  mirror. 

The  superposition  principle  can  be  used  to  describe  reflection 
from  more  complex  mirror  surfaces.  In  case  of  a  mirror  with  a  hole  the 
surface  integral  (7)  over  S  is  expressed  as  /g  =  fg^  -  where  is 
defined  by  the  mirror  surface  including  the  hole  surface,  and  Sj  is  the 
surface  of  the  hole  only.  The  total  reflection  matrix  R  is  given  by  R  = 
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R(Sj)  -  R(S7),  the  difference  in  the  reflection  matrices  associated  with 
mirrors  and  S2  respectively.  The  transmission  matrix  T  through  a 
screen  with  an  aperture  of  area  S  is  given  by  T  =  -  R,  R  being  the 
reflection  matrix  for  a  mirror  matching  the  aperture  S.  The 
transmission  matrix  for  radiation  diffracted  behind  a  finite  size  mirror 
is  given  by  T'  =  1  -  elrt  R  where  1  is  the  identity  matrix. 


V.  REFLECTION  OF  THE  LOWEST  ORDER  MODE 

The  computation  of  all  the  truncated  integrals  for  finite  mirror 

surface  is  nontrivial.  Most  applications  however  involve  the  (0,0) 

lowest  order  mode  as  the  dominant  mode  in  both  incoming  and  reflected 

radiation.  The  strategy  here  is  to  compute  the  element  Rqq  of  the 

reflection  matrix  first.  Then  the  waist  for  the  reflected  modes  W  can 

o 

be  selected  so  that  it  maximizes  Rqq.  The  optimum  representation 
condition 


9R, 


00 

00 


9a 


=  0, 


(26) 


puts  the  maximum  amount  of  the  reflected  radiation  in  the  lowest  mode  (a 
different  mode  and  matrix  element  may  be  chosen,  if  desired).  It  is 
pointed  out  that  (26)  does  not  improve  the  properties  of  the  reflected 
radiation.  It  enables  one  to  choose  the  best  representation  in  terms  of 
minimizing  the  coefficients  of  the  undesired  modes  for  the  scattered 
radiation.  Once  Wq  is  fixed  by  (26)  then  the  exact  location  and  size  of 
the  waist(s)  for  the  reflected  modes  is  determined  by  solving  the  system  of 
equations 


W 


0 


w 


o 


1/2 

f 


1 

R 


0 


l2+b2 
o  o 


(27) 


» 


14 


The  element  is  identical  in  both  representations  since  the  lowest 
order  mode  Uqq  is  the  same  in  rectangular  and  cylindrical  coordinates. 
Performing  the  integration  (20)  yields  R^q  to  first  order  in  £ 


,00 

*00 


2  a 


1+a 


1  -  e 


/i  2.  2 

-(1+a  )y 


4  oa2). 


(28) 


Note  that  the  first  order  term  vanishes  and  the  lowest  correction  is  of 

2 

second  order  in  5  •  The  exact  dependence  on  the  mirror  size  p  is 
parametrized  by  y  =  p/Vo>  while  a  =  Wj/Vq  parametrizes  the  ratio  of  the 
incoming  and  scattered  radiation  spot  sizes  at  the  mirror.  The  optimization 
condition  3R°°(0)/3a  =  0  yields,  a2  =  1  +  exp[-(l+a2)y2  ]  [a4+(2y2+l)a2-l) . 
In  case  the  mirror  cross  section  is  much  larger  than  the  spot  size  of  the 
incoming  mode,  y  »  1,  a  -*  1  and  the  reflected  spot  size  at  the  mirror 
matches  the  incoming,  Wq  *  V^. 

Large  mirror  size  is  desired  to  maximize  the  total  reflection 

coefficient.  The  reflection  coefficient  0^  is  given  by  =  P</Pi  wllere  t*le 

i  2  i  2 

incoming  flux  is  P.  =  |c  |  =  E  |c  |  and  the  outgoing  flux  is  given  by 


Po=  k 


o  i2 


«  •  *'l2  ■  Z  E  I"™  c*Ml: 

mn  pq 


(29) 


In  Fig.  2  we  plot  for  the  lowest  order  incoming  mode  as  a  function  of 
y'  =  cos($/2)  p/WQ  =  cos($/2)  y.  y'  parametrizes  the  size  of  the  mirror 
surface  projection  into  the  plane  perpendicular  to  the  incoming 

radiation  direction.  The  incoming  radiation  has  a  wavelength  X  =  ly 

-4  -2  3 

(10  cm),  waist  w^  =  10  cm  at  distance  1^  =  2.5x10  cm  from  the  mirror 

3 

and  radius  of  curvature  (at  the  mirror)  Rj  =  2.5x10  cm.  The  mirror  has 

3 

a  radius  of  curvature  R  =  2.5x10  cm,  yielding  reflected  modes  of  R  = 

in  o 


15 


3 

2.5x10  (again  1  and  wq  depend  on  the  choice  of  VQ).  In  Fig.  3  we  plot 
the  reflection  coefficients  of  the  lowest  order  mode  (0,0)  into  the 
first  25  modes  (p,q)  with  p  <  q  <  5,  as  a  function  of  p' .  The 
deflection  angle  is  90°  and  the  ratio  of  the  spot  sizes  is  1. 
Increasing  mirror  size  maximizes  the  diagonal  element  and  minimizes 
scattering  into  other  modes. 

In  Fig.  4  we  fix  the  mirror  size  p'  =  2  and  the  angle  ♦  =  90°  and 

vary  the  spot  size  ratio  a.  The  best  representation,  maximizing 

and  minimizing  R^®  is  obtained  at  a  =  1.  However,  for  small  mirror  p'  = 

pq 

0.66,  the  maximum  for  Rqq  occurs  at  a  «  0.70  (see  Fig.  5).  Radiation 
reflected  off  mirrors  smaller  than  the  incoming  spot  size  is  best 
described  by  outgoing  modes  of  reduced  spot  size  Wq  <  V. .  Also  note 
from  Fig.  5b  that  for  small  mirror  size  the  total  power  reflected  into 
the  first  25  modes  never  exceeds  802  of  the  incoming  flux;  even  with 
many  more  modes  remains  less  than  1.  In  Fig.  6  the  reflection 
coefficients  R®^  are  plotted  as  functions  of  the  angle  of  deflection  ♦ 

pq 

for  fixed  a  =  1,  p'  =  2.  It  is  seen  that,  for  sufficiently  large 

reflecting  surface  and  good  choice  of  the  spot  size  Vo,  the  reflection 
matrix  is  not  very  sensitive  to  $.  Comparison  of  Figs.  2  and  3  with  the 
rest  of  the  plots  shows  that  the  relative  mirror  size  to  the  radiation 
spot  size  is  the  most  important  parameter  to  determine  the  reflection 
into  other  than  the  incoming  modes. 
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Figure  Captions 


Figure 

Figure 

Figure 


Figure 

Figure 

Figure 


Reflection  geometry. 

Plot  of  the  total  reflection  coefficient  for  the  lowest  order 
mode  as  a  function  of  the  mirror  size  y'  for  $  =  90°  and  Wq  =  V^. 
Plots  of  the  reflection  matrix  elements  of  the  lowest  order  mode 
(0,0)  into  the  first  25  modes  p  <  q  <  5.  The  magnitude  |R^j  |  is 
plotted  against  the  relative  mirror  size  y' .  The  angle  of 
deflection  $  =  90°  and  a  =  1  (W^  =  VQ). 

Plots  of  (a)  the  reflection  matrix  elements  l^pql'  an(*  (b)  the 
reflection  coefficient  it^  into  the  first  25  modes,  against  the 
spot  size  ratio  a  for  y'  =  2  and  $  =  90°. 

Same  as  in  Fig.  4  for  y'  =0.5. 

Plots  of  the  reflection  matrix  elements  lRpq|  against  the  angle  of 
deflection  ♦  for  y'  =  2  and  a  =  1. 
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